Where are the black hole entropy degrees of freedom ? 
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Understanding the area-proportionality of black hole entropy (the 'Area Law') from an underlying 
fundamental theory has been one of the goals of all models of quantum gravity. A key question 
that one asks is: where are the degrees of freedom giving rise to black hole entropy located? Taking 
the point of view that entanglement between field degrees of freedom inside and outside the horizon 
can be a source of this entropy, we show that when the field is in its ground state, the degrees of 
freedom near the horizon contribute most to the entropy, and the area law is obeyed. However, 
when it is in an excited state, degrees of freedom far from the horizon contribute more significantly, 
and deviations from the area law are observed. In other words, we demonstrate that horizon degrees 
of freedom are responsible for the area law. 



PACS numbers: 04.70.Dy, 03.67.Mn, 03.65.Ud,05.50.+q 

The area proportionality of black hole entropy (the 
'Area Law' (AL)), 



(l p = y/hG/c 3 = Planck length) (1) 



which differs from the volume proportionality of famil- 
iar thermodynamic systems, has been conjectured to be 
more fundamental in some sense (the Holographic Hy- 
pothesis) . Black holes are also regarded as theoretical lab- 
oratories for quantum gravity. Thus, candidate models of 
quantum gravity, such as string theory and loop quantum 
gravity, have attempted to derive the 'macroscopic' AL 
(1) by counting 'microscopic' degrees of freedom (DOF), 
using the Von-Neumann/Boltzmann formula [1]: 



S = -Tr[p\n{p)] = InO, 



k B = 1 



(2) 



where p and n correspond to the density matrix and num- 
ber of accessible micro-states, respectively. Depending on 
the approach, one either counts certain DOF on the hori- 
zon, or abstract DOF related to the black hole, and there 
does not appear to be a consensus about which DOF are 
relevant or about their precise location [2]. 

In this article, we attempt to answer these questions 
in a more general setting, which in fact, may be relevant 
in any approach. We adopt the point of view that the 
entanglement between quantum field DOF lying inside 
and outside of the horizon leads to black-hole entropy. It 
was shown in Refs. [3, 4] (see also, Refs. [5, 6]) that the 
entanglement entropy of a massless scalar field propagat- 
ing in flat space-time (by tracing over a spherical region 
of radius R) is proportional to the area of the sphere 
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is the UV cutoff . 



(3) 



This suggests that the area law is a generic feature of en- 
tanglement, and acquires a physical meaning in the case 
of black-holes, the latter's horizon providing a natural 
'boundary' of the region to trace over. Note that Eq.(3) 



will continue to hold if the region outside the sphere is 
traced over instead. 

The relevance of S oat to S m can also be understood 
from the fact that both entropies are associated with the 
existence of the horizons [7]. Consider a scalar field on 
a background of a collapsing star. At early times, there 
is no horizon, and both the entropies are zero. How- 
ever, once the horizon forms, S BH is non-zero, and if the 
scalar field DOF inside the horizon are traced over, S f ont 
obtained from the reduced density matrix is non-zero as 
well. 

In Refs. [3, 4], along with the fact that calculations 
were done in flat space-time, a crucial assumption was 
made that the scalar field is in the Ground State (GS). In 
Refs. [8, 9], the current authors studied the robustness of 
the AL by relaxing the second assumption, and showed 
that for Generic Coherent States (GCS) and a class of 
Squeezed States (SS), the law continues to hold, whereas 
for the Excited States (ES), one obtains: 



R 



0(1) 



(4) 



where a < 1, and higher the excitation, the smaller its 
value. In this article, we attempt to provide a physical 
understanding of this deviation from the AL, by show- 
ing that for ES, DOF far from the horizon contribute 
more significantly than for the GS. We also extend our 
flat space-time analyses for any (3 + 1)— D spherically 
symmetric non-degenerate black-hole space-times [10]. 

We begin by considering a scalar field 4>{x) propagating 
in a Schwarzschild space-time [17]: 



■f(r)dt 2 + j^+r 2 dQl 



m = i- r i 

r 



(5) 



where r is the horizon. Transforming to Lemaitre coor- 
dinates (t,R) [12]: 



t + r Q 
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the line-element (5) becomes: 
~3{R-t 



ds 2 — — dr 2 + 



ro 



dR 2 + 



3{R 
2 



ro 



r^dQ 2 . (7) 



Note that R(t) is everywhere space- like (time- like), the 
metric is non-singular at r = ro (corresponding to 3(i? — 
r)/2ro = 1) and the metric is explicitly time-dependent. 
The Hamiltonian of a massless, minimally coupled scalar 
field in the background (7) is given by 



lm Jr 



oo 

dR 



2P l j n (r,R) 3r 
3(R-r) + Y {R - T) 



(8) 



x [d R( p lm (r,R)} 2 + ? -H-l(t+l)tf m {T,R) 



where P; m (r, R) is the canonical conjugate momenta of 
the spherically reduced scalar field <pi m (r, R) , such that 
[4>i m {R, To),Pi m (R', t )] = i5(R—R') and I denotes angu- 
lar momenta. Note that, in these coordinates, the scalar 
field and its Hamiltonian explicitly depend on time. 

Next, choosing a fixed Lemaitre time (say, t = 0) and 
performing the following canonical transformation [16] 

Pi m (r)=Vr~Tr lm (r) , 0, m ( r ) = , (9) 

the Hamiltonian (8) transforms to: 



H(0) 



re 



dr 
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irim{r)+r [d r —\ + ^ fim(r) 

(10) 

which is simply the Hamiltonian of a free scalar field in 
flat space-time! Eq.(10) holds for any fixed t, for which 
the results of Refs. [4, 9], namely relations (3) and (4), 
go through, provided one traces over either the region 
R £ [0, |r ) or the region R £ [|r ,oo) [10]. Extension 
to any non-degenerate spherically symmetric space-times 
is straightforward [18]. 

Remaining steps in the computation of 5 ent are: 
1. Discretize the Hamiltonian (10), i. e., 



1 N 



IT 2 

Im.j 



1 \ 2 / fhn.J flm,j + l ^ 



+ -T2 



(11) 



where iri m ,j denotes the conjugate momenta of ipi m ,j , L = 
(N + l)a is the box size and a is the radial lattice size. 
This is of the form of the Hamiltonian of N coupled HOs 



N N 

ff =oE^+9 E x * K U x i i,3 = l,---,N (12) 



with the interaction matrix elements given by: 



KlJ ~ i 2 
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1(1 + 1) Sij + - SnSji + [ N 




i,j(i^l,N) 



(13) 



where the last two terms are the nearest-neighbor inter- 
actions. 

2. Choose the state of the quantum field: For GS (r = 
l,ai = 0), GCS (r = l,a>i = arbitrary), or a class of 
SS (a — 0, r = arbitrary), the N-particle wave function 
ip(xi . . . xn) is given by [9] 



ip(xi ■ --x N ) = 










exp 



1 1/2 , n2 



.(14) 



For ES (or 1-Particle state), the N-particle wave function 
i\)(x\ . . . Xn) is given by 



ip{x\ ■ ■ ■ x N ) 



2^ 1 1 v- A 

! i — 1 



4E4 



(15) 



i=l 



where a T = (ai, . . . , ajv) are the expansion coefficients 
(normalization requires a T a = 1). For our computa- 
tions, wc choose a T = l/^/o(0, • • • 0, 1 • • • 1) with the last 
o columns being non-zero. For details, see Ref. [9]. 
3. Obtain the density matrix: For an arbitrary wavefunc- 
tion ip(xi, . . . ,Xn), the density matrix — tracing over 
first n of the N field points — is given by: 

/n 
( [ dxi ip(xi,...,x n ;ti,...,t N - n ) 
»=i 

x ^*(o;i,...,a; n ;ti,...,^_ n ) (16) 

where: tj = x n+j ,t = ti, . . .,t N - n , 3 = l-(N -n),x T = 
(xi, ...,x n ;ti,... ,t N -i) = (xi, . . .,x n ;t). This step, in 
general, cannot be evaluated analytically and requires 
numerical techniques. For GS/CS/SS, substituting (14) 
in the above expression and using the relation (2) gives 
Eq. (3). For ES, on the other hand, this leads to the 
relation (4) [9]. 

Now, to locate those DOF which are responsible for 
entropy, we take a closer look at the interaction matrix 
(13). As mentioned earlier, the last two terms are the 
nearest-neighbor interactions between the oscillators and 
constitute entanglement. As expected, if these terms are 
set to zero (by hand), S^nt vanishes. Instead, let us do 
the following: 

(i) Set the interactions to zero (by hand) everywhere ex- 
cept in the 'window' (q— s < i < q+s), with s < q < n—s, 
i.e., restrict the thickness of the interaction region to 
2s + 1 radial lattice points, while moving it rigidly across 
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from the origin to horizon. In Fig. (1), we have plot- 
ted the percentage contribution to the total entropy as a 
function of q for the GS/GCS/SS (o = 0, solid thin curve) 
and the ES (o = 30(50), bold (light) thick curve). [We 
choose N = 300 and n = 100, 150, 200. The numerical 
error in the evaluation of the entropy is less than 0.1%.] 
Fig. (1) shows that: 

• When the interaction region does not include the 
horizon, the entanglement entropy is zero and it is 
maximum if the horizon lies exactly in its middle. 
The first observation confirms that entanglement 
between DOF inside and outside the horizon con- 
tributes to entropy, while the second suggests that 
DOF close to the horizon contribute more to the 
entropy compared to those far from it. 

• When the number of excited states is increased (i. 
e. o = 30,50), the percentage contribution to the 
total entropy close to the horizon is less compared 
to that of GS and the (bold/light thick) curves are 
flatter. These indicate that, for ES, there is a sig- 
nificant contribution from the DOF far-away from 
the horizon. 
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entropies [S en t(t)} vs t for n — 100, 150 and 200. [Here 
again, the solid thin curve is for GS and bold (light) thick 
curve for o = 30(50).] We infer the following: 

• For the GS/GCS/SS, the entropy reaches the GS 
entropy for as little as t = 3. In other words, the 
interaction region encompassing DOF close to the 
horizon contribute to most of the entropy for the 
GS/GCS/SS. 

• In the case of ES, the entropy reaches the ES en- 
tropy only for t = 15(20) [for o = 30(50)]. This in- 
dicates that: (a) The DOF far-away from the hori- 
zon have a greater contribution than that of GS and 
(b) As the number of excited states increases, con- 
tribution far-away from the horizon also increases. 

This leads us to one of the main conclusions of 
this article: the greater the deviation from the AL, 
the larger is the contribution of the DOF far-away 
from the horizon. It can be shown that the density 
matrix for the ES is more spread out than for the 
GS. 
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FIG. 1: Plot of the percentage contribution to the entropy for 
the GS and ES as we move the window q — 2 < i < q + 2 from 
q = 2 to q = n, for N = 300 and n = 100, 150 and 200. The 
solid thin curve is for GS o = and bold (light) thick curve 
for o = 30(50). 

(ii) To further investigate, we now set the interactions 
to zero (by hand) everywhere except in the window p < 
i < n, with the horizon as its outer boundary, and vary 
the width of the window t = n — p from to n. For 
t = n, the total entropy is recovered, while for t = 0, 

1. e. zero interaction region, entropy vanishes. In Fig. 

2, we have plotted the normalized GS/GCS/SS and ES 



FIG. 2: Plot of the percentage contribution of Sent(t) for the 
GS and ES. Here again, N = 300, n = 100, 150 and 200, and 
the solid thin curve is for GS and bold (light) thick curve for 
o= 30(50). 

(iii) To understand this further, let us define 

Sent (t) 



Apc(t) = pc(t) — pc(t — 1) where pc(t) 



S c 



x 100. 



where pc(t) is the percentage contribution to the total en- 
tropy by an interaction region of thickness t and Apc(t) 
is the percentage increase in entropy when the interac- 
tion region is incremented by one radial lattice point. In 
other words, Apc(t) is the slope of the curves in Fig. 
(2). In Fig. (3), we have plotted Apc(t) vs (n — t) for 
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GS and ES. For the GS/GCS/SS, it is seen that when 
the first lattice point just inside the horizon is included 
in the interaction region, the entropy increases from 
to 85% of the GS entropy. Inclusion of the next three 
points add another 9%, 3%, 1% respectively. The contri- 
butions to the entropy decrease rapidly and by the time 
the (n/3) th point is included, entropy barely increases by 
a hundredth of a percent! For ES however, inclusion of 
one lattice point adds 70(50)%, for o = 30(50), to the 
entropy, while the next few points add 9%,4%,3%--- 
respectively. The corresponding slopes are smaller. 
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FIG. 3: Plot of Apc(t) vsn-t for GS and ES. Here again, 
N = 300, n = 100, 150 and 200, and the solid thin curve is 
for GS and bold (light) thick curve for o = 30(50). 

Figures (2) and (3) suggest the next key result of this 
article - most of the entropy comes from the DOF close 
to the horizon. However, the DOF deep inside must also 
be taken into account for the AL (3) to emerge for the 
GS/GCS/SS, and the law (4) to emerge for the ES. These 
DOF affect the horizon DOF via the nearest-neighbor 
interactions in (13). 

Our work clearly demonstrates the close relationship 
between the AL and the horizon DOF, and that when 
the latter become less important, the entropy scales as a 
power of area less than unity. This can be understood 
from the following heuristic picture: taking the point of 
view that S BH is proportional to A/", the number of DOF 
on the horizon, and that for the GS/CS/SS there is one 
DOF per Planck-area, such that Af cx Ah, the AL follows 
(This is known as the it-from-bit picture [14]). For ES 
however, since this number is seen to get diminished, it 
will be given by another function Af = f(A) < A. Now, 
since current results are expected to be valid when i> 1 
(in Planck units) (such that backreaction effects on the 
background can be neglected), it is quite plausible that 
f(A) ~ A a , a < 1, just as we obtain. Note that the 



above reasoning continues to hold when the outside of the 
horizon is traced over. Another way of understanding our 
result is as follows: all interactions being of the nearest- 
neighbor type [cf. Eq.(13)], the degrees of freedom deep 
inside the horizon influence the ones near the horizon 
(and hence contribute to the entropy) only indirectly, i.e. 
via the intermediate degrees of freedom. Evidently, their 
influence diminishes with their increasing distance from 
the horizon. When they are excited however, they have 
more energy to spare, resulting in an increased overall 
effect. Further investigations with superpositions of GS 
and ES are expected to shed more light on this [15]. 
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